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ABSTRACT 

Gas motions in the hot intracluster medium of galaxy clusters have an important 
effect on the mass determination of the clusters through X-ray observations. The cor- 
responding dynamical pressure has to be accounted for in addition to the hydrostatic 
pressure support to achieve a precise mass measurement. An analysis of the velocity 
structure of the ICM for simulated cluster-size haloes, especially focusing on rotational 
patterns, has been performed, demonstrating them to be an intermittent phenomenon, 
strongly related to the internal dynamics of substructures. We find that the expected 
build-up of rotation due to mass assembly gets easily destroyed by passages of gas- 
rich substructures close to the central region. Though, if a typical rotation pattern is 
established, the corresponding mass contribution is estimated to be up to ~ 17% of 
the total mass in the innermost region, and one has to account for it. Extending the 
analysis to a larger sample of simulated haloes we statistically observe that (i) the dis- 
tribution of the rotational component of the gas velocity in the innermost region has 
typical values of ~ 200 — 300km/s; (ii) except for few outliers, there is no monotonic 
increase of the rotational velocity with decreasing redshift, as we would expect from 
approaching a relaxed configuration. Therefore, the hypothesis that the build-up of 
rotation is strongly influenced by internal dynamics is confirmed, and minor events 
like gas-rich substructures passing close to the equatorial plane can easily destroy any 
ordered rotational pattern. 
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1 INTRODUCTION 

Within the hierarchical structure-formation scenario, galaxy 
clusters are key targets that allow us to study both the dy- 
namics on the gravity-dominated scale and the complexity 
of astrophysical processes dominating on the small scale. In 
such studies their mass is one of the most crucial quanti- 
ties to be evaluated, and the bulk properties measured from 
X-ray observations still provide the best way to estimate 
the mas s, primarily on the assumption of hydrostatic equi- 
librium dSarazinlllQSi l. Mass estimates rely then on the as- 
sumptions made about the cluster dynamical state, since 
the Hydrostatic Equilibrium Hypothesis (HEH) implies that 
only the thermal pressure of the hot ICM is taken into 
account asia. Tormen fc Moscardinil 120041 '). Lately, it has 
been claimed in particular that non-thermal motions, as ro- 
tation, could play a significa nt role in supporting the ICM 
in the innermost region (e.g. iLau. Kravtsov fc Nagaill2009l : 
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iFang. Humphrey fc Bu"ot3l2009D biasing the mass measure- 
ments based on the HEH. The analysis of simulated cluster- 
like objects provides a promising approach to get a better 
understanding of the intrinsic structure of galaxy clusters 
and the role of gas dynamics, which can be eventually com- 
pared to X-ray observations. Because of the improvement 
of numerical simulations, the possibility to study in de- 
tail the physics of cluster s has enormously increased (see 
Borgani fc Kravtsov 2009, for a recent comprehensive re- 
view) and future satellites dedicated to high-precision X- 
ray spectroscopy, such as ASTRO-H and IXO, will allow 
to detect these ordered motions of the ICM. With this per- 
spective, we perform a preliminary study on the ICM struc- 
ture for some clusters extracted from a large cosmological 
hydrodynamical simulation, investigating in particular the 
presence of rotational motion in the ICM velocity field. 

The paper is organized as follows. We describe the nu- 
merical simulations from which the samples of cluster-like 
haloes have been selected in Section (2] In Section [3] we con- 
sider a first set of simulated clusters and present results on 
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build-up of rotation in the halo core for single cases of study 
(Section 13.11 and Section 13. 2p , while results about the con- 
tribution to the mass estimations are given in Section [H A 
second sample of clusters is then statistically investigated in 
Section [5] We discuss our results and conclude in Section |6l 
Appendix A is devoted to comment on the effects of ar- 
tificial viscosity, while in Appendix B we briefely comment 
on the ellipticity profiles of the simulated clusters. 



2 NUMERICAL SIMULATIONS 

We consider two sets of cluster-like haloes selected from two 
different parent cosmological boxes. In both cases the cosmo- 
logical simulations we re performed with the TreePM/SPH 
code GADG ET-2 ijSpringel. White fc HernauistI 120011 : 
ISpringell [20051 ) ■ assuming a slightly different cosmological 
model (a standard ACDM model and a WMAP3 cosmology, 
respectively) but including t he same physical processe s 
governing the gas component (|Springel fc Ifernauistll2"003h . 
i.e. radiative cooling, st ar formation, and supernova feed- 
back {csf simulation, see lDolag et al.l (|2009 ) and references 
therein for a more detailed overview on different runs of 
the parent hydrodynamical simulations we refer to in our 
work). Additionally, we refer to simulations of the same 
objects without including radiative processes as ovisc. 

Set 1. The first data set considered consists of 9 
cluster-size haloes, re-simulated with higher resolution 
using the "zoomed initi al con dition" (ZIC) technique 
ijTormen. Bouchet fc White! 1 19971 ). The clusters have 
been originally extracted from a large-size cosmological 
simulation of a ACDM universe with ilo ~ 0.3, h = 0.7, 
(J8 = 0.9 and rib = 0.039, within a box of 479/i"^Mpc 
a side. The final mass-resolution of these simulations is 
moM = 1.13 X lO^/i~^M0 and m,gas = 1.69 x lO^/i'^M©, 
for the DM and gas particles, respectively. The spatial 
resolution for Set 1 reaches 5/i~^kpc in the central parts 
and for the most massive clusters we typically resolve up 
to 1000 self-bound sub-structures within R^r (as shown 
m iDolag et al.|[2009l ). The main haloes have masses larger 
than ~ 1.1 X 10^'*/i~^Mq and they have all been selected in 
a way that they are quite well-behaved spherically-shaped 
objects at present epoch, although a fair range from isolated 
and potentially relaxed objects to more disturbed systems 
embedded within larger structures is available. Having a 
reasonable dense sample within the time domain (e.g. 50 
outputs between z = \ and today) and high resolution, 
we can study how common and significant the rotational 
support of the ICM is. In particular we focus on the detailed 
evolution of such rotational motions. 

Set 2. In the second data set we analyzed a volume 
limited sample of cluster-size haloes, where we computed 
the distribution of the ICM rotational velocity and com- 
pared their distribution at different redshift. This second set 
of clusters has been extracted from a large size cosmological 
simulation with a box-size of 300/i~^Mpc, simulated with 
2 X 768^ particles, assuming the 3-year WM AP values 
for the cosmological parameters (Spergel et alJ 20^), i.e. 
fio = 0.268, S^Jfc = 0.044, erg = 0.776 and h = 0.704. The 
final mass-resolution for this second set of simulations is 
niBM = 3.71 X IO^/i^^Mq and rugas = 7.28 x lO*/i"^M0. 
Given the larger sample, a fair investigation of the amount 
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Figure 2. Main halo mass accretion history for the four most 
massive clusters in Set 1: g51 {solid, black line), gS {dotted, blue), 
gl {dashed, green), and g72 {dot-dashed, light-blue). The mass 
of the main halo is plotted in units of lO^'*/i~^M0 as function of 
redshift. 

of rotational support within the ICM from a statistical 
point of view is then possible. 

The two sets of simulations analyzed were performed 
with a different value of as, meaning that differences in the 
merging histories can be introduced. Nevertheless, we stress 
that the purpose of the second set is only to enlarge the 
statistics on the build-up of rotational motions in simulated 
galaxy clusters, and not a direct comparison of single 
objects to the objects of Set 1. Therefore, we are confident 
that our conclusions do not depend on these differences in 
the parameters of the two simulations. 



3 VELOCITY STRUCTURE OF THE ICM 

In the first place, we have been investigating the velocity 
structure of the ICM for all the main haloes of the Set 1, 
looking for evidence of rotational patterns in the most cen- 
tral regions of the simulated clusters, especially for the most 
massive haloes. 

Although a variety of objects is offered (ranging from 
more complex structures sitting in a denser environment to 
quite isolated haloes), all the cluster- like haloes have been 
originally selected to be fairly regularly shaped. For this as- 
pect, the whole Set 1 is biased towards quite relaxed clus- 
ters and we would expect to find a more significant amount 
of rotation than on average. Among the four most massive 
systems, for which Fig. [T] gives the X-ray surface bright- 
ness maps in the three projected directions, we identify as 
the most relaxed clusters g51 and g8, while gl and g72 are 
disturbed systems still suffering at the present epoch from 
recent major mergers 

An indication of this differences is clearly shown in 
Fig. [21 where we plot the mass accretion history for the 
main progenitor of each cluster as function of redshift. The 
curves referring to g51 and g8 (solid black line and dotted 
blue line, respectively) show a smoother mass assembly 
(at late times, i.e. z < 0.3) if compared with those for gl 
and g72 (dashed green line and dot-dashed light-blue line. 
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Figure 1. X-ray surface brightness maps along tlie tliree projection axis for the four most massive cluster-like haloes in Set 1 (from top 
to bottom: g51, g8, gl, g72). Each map is 2Mpc-side, enclosing therefore the region of about ~ Rsoo- 
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respectively), whose curves show bumps related to major 
merging events down to very low redshift. In addition, a 
visual inspection of the X-ray surface brightness maps 
presented in Fig. [1] definitely suggests g51 to be the less 
sub-structured halo. 

In this perspective, we address g51 as the best case of 
study to explore the build-up of rotational motions in the 
cluster central region as a consequence of the cooling of the 
core. Among the other clusters also gl shows interesting 
features in its velocity field that are worth to be investigated 
in more detail and compared to the case of g51 to better 
characterize the occurrence of ordered rotational motions 
in the intracluster gas (see Appendix |B] for a comment on 
the isophote ellipticities of g51 and gl among the four most 
massive haloes from Set 1 and how the gas shapes relate to 
intrinsic rotational gas motions). 



plane" on which we can calculate the tangential component 
of the velocity, vtan- This plane is taken to be perpendicu- 
lar to the direction of the mean gas angular momentum, j, 
calculated averaging over the gas particles within the region 
where we want to investigate the rotational motion within 
the ICM, i.e. < O.liisoo (here, the overdensity of 500 is de- 
fined with respect to the critical density of the Universe). 
This definition of the "best equatorial plane" allows us to 
emphasize and characterize the rotation of the gas in an ob- 
jective way for all clusters, whenever it appears. 

To perform our analysis, we rotate the halo such that 
the new z-axis is aligned with the direction of j and the 
new xy plane easily defines the best equatorial plane. After 
subtracting an average bulk velocity for the gas component 
within the region corresponding to -R500, we compute the 
tangential component of the velocity on this plane. We con- 
sider a 40/i~^kpc slice of the simulation box containing this 
plane for all the calculations hereafter. 



3.1 Rotational patterns in the ICM 

The two cases analyzed in detail (the isolated regular clus- 
ter, g51, and the disturbed massive halo gl) are particularly 
interesting for our purpose, since their velocity structure at 
redshift z = shows two opposite pictures, namely strong 
rotational patterns for gl and almost no gas rotation for 
g51. 

In the classic cooling flow model, gas rotation is ex- 
pected near the center of the flow because of mass and 
angu lar momentum conservation (e.g. Mathews & Brighenti 
I2OO3I ). Though, the rate of cooling gas predicted by the clas- 
sical paradigm of cooling fiows is rarely observed in real clus- 
ters, implying that feedback processes must p lay a role in 
preventing cooling (jMcNamara fc NulsenlfioOTI ) . In contrast, 
in simulations, the strong cooling in the central region of the 
simulated cluster-like haloes suggests that relaxed objects 
should build up significant rotational motion in the inner- 
most region where gas is infalling and contracting under the 
conservati on of angular momentum. A s reported in the liter- 
ature (e.g. iFang. Humphrev fc Buo t3 i2009.1 this effect is ex- 
pected to be particularly evident in simulated clusters that 
can be identified as relaxed objects. Therefore it is interest- 
ing to point out that, in our sample, not even in the object 
with the smoothest accretion history and less substructured 
morphology significant rotational patterns establish in the 
ICM velocity structure as a consequence of collapse. 

Dealing with hydrodynamical simulations though, the 
build-up of rotation in the central region of clusters can 
be also enhanced by an excess of gas cooling that has 
been found to overprod uce the observed cosmic abundance 
of st ellar material (e.g iKatz fc Whiti 1 19931 : iBalogh et al] 
[200I) in absence of very strong, not yet fully under- 
stood feedback processes. In our simulations, the imple- 
mentation of a multi-phase model fo r star formation (e.g. 
Katz. Weinberg fc Hernguistl Il996l : ISpringel fc Hernguistl 
20031 ) and the treatment of the thermal feedback pro- 
cess, including also galactic winds associated to star for- 
mation, is able to pa rtially reduce the over-cooling problem 
l|Borgani et al.|[2006l ) . This fact plausibly contributes to pre- 
vent significant rotation. 

In order to study in detail the rotational component of 
the ICM velocity for a halo, we first define a "best equatorial 



3.2 A case study: g51 vs. gl 

As a case study, we particularly focus on g51, an iso- 
lated massive cluster with gravitational mass of M2oom ~ 
1.34 X 10^^/i"^Mq and a size of J?200m = 2.28ft"^Mpc at 
z = 0, and we compare it with the other extreme case men- 
tioned, gl, which is instead a strongly disturbed system with 
Af2oom = 1.49 X lO^^/i"^M0 and R2oom = 2.36/i"^Mpc. 
Here, -R200m is defined as the radius enclosing the region 
with density equal to 200 times the mean density of the Uni- 
verse, and M2oom is the mass within i?200m ■ M2oom has been 
used throughout our study as reference quantity to select 
haloes, but we always carry out our calculations by referring 
to iisoo and M500, where the overdensity of 500 is instead 
defined with respect to the critical density of Universe, mo- 
tivated by a possible comparison to real X-ray observations. 
In the case of g51 and gl, we have i?5oo = 1.09, 1.20/i~^Mpc 
and M500 = 7.46 x 10^"*, 9.98 x 10"/i"^MQ respectively, at 
z = 0. 

From the considerations made in Section |3] about its 
shape and accretion history, g51 is likely to be, in a global 
sense, the most relaxed object in the sample. In spite of 
this, at 2 = the velocity structure of the ICM in the inner- 
most region is far from showing a clear rotational pattern 
as expected from a nearly homogeneous collapse process. 
However it shows some rotational pattern at intermediate 
redshift. 

In Fig.|3]we plot the rotation velocity profile, Vtan{r), for 
the two interesting cases (at z = 0) out to -Rsoo- As explained 
in the previous Section, vtan is the tangential component of 
the ICM velocity, calculated in the best equatorial plane. In 
order to compute the radial profile displayed in Fig. O we 
make use of radial bins in the plane to calculate the mass- 
weighted average value of vtan of the gas particles at each 
r. We have chosen 14/i~^kpc, as optimal bin width on the 
base of both resolution and statistical motivation. 

The radial profiles for the rotational component of the 
gas velocity reflect the presence of a non-negligible rota- 
tional pattern in the ICM of the disturbed system, while no 
significant rotation is built up in the relaxed one. The profile 
of g51 (left panel in the Figure, solid curve) shows relatively 
low values at small radii, and increases significantly only at 
radii larger than 0.3i?5oo, where the rotational compo- 
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Figure 3. Rotational velocity as function of the radius out to -R500 for a relaxed cluster (g51, left panel) and for a highly disturbed 
system (gl, right panel), at 2; = 0. The tangential component of the ICM velocity is calculated on the best equatorial plane, i.e. the 
plane perpendicular to the direction of the mean gas angular momentum in the region within O.liJsoo- Two runs are compared: the csf 
simulation {solid line), including radiative cooling, star formation and supernova feedback, and the ovisc simulation {dashed line), in 
which all these physical processes are omitted and only a treatment for artificial viscosity is considered. 



nent of the velocity is likely to be dominated by some bulk 
rotational motions, plausibly related to a subhalo orbiting 
in the main halo close to Rsoo- The value of vtan decreases 
instead with increasing r for gl (right panel in the Figure, 
solid curve), where the rotational velocity reaches almost 
lOOOkm/s in the innermost region. Also, it is interesting to 
compare with the rotational velocity profile for two counter- 
part haloes, simulated without including star formation and 
cooling (dashed curv es). In such simul ation, referred to as 
ovtsc simulation (see lDolag at alJl2009l . as an overview) the 
overcooling problem is completely avoided because no stars 
are formed at all, and no significant rotation is expected to 
build up in the center of the cluster-like haloes. Let us note 
that for g51, the curves referring to the two simulations have 
a significantly similar trend, while for gl the csf simulation 
(solid curve) and the non-radiative one diverge towards the 
center, increasing in the former and decreasing in the latter. 
While major events occurring close to « = in the merging 
history of gl could explain the high values found for vtan 
in the innermost region, no major mergers happen to char- 
acterize the history of g51 at late time. Therefore a further 
zoom onto g51 is required in order to understand the de- 
tails of the processes that lead to the build-up or to the 
disruption of gas rotation in the halo core. 

Rotational velocity evolution. The possibility to 
track back the history of the cluster-size haloes given by 
simulated data, allows us to follow the redshift evolution 
of the rotational component of the ICM velocity in the 
innermost region of g51, taken to be O.lTisoo. Up to z — 2, 
a mass-averaged value of the tangential component of the 
ICM velocity has been calculated in the best equatorial 
plane, so that rotation can be emphasized best whenever 
there is one. At each redshift, the orientation of the best 
equatorial plane has been adjusted to be perpendicular 
to the direction of the mean gas angular momentum, as 
previously defined. 



While in the literature w e find an inspiring work 
(e.g. iFang. Humphrev fc Buotj I2OO9I ) where values for 
the rotational velocity in the central region of a re- 
laxed cluster-like halo rise above lOOOkm/s, in our 
study this never happens and values generally increase 
up to 650km/s as a maximum, except for high peaks 
probably related to major merging events. These differ- 
ences are likely to be related to the different amount of 
baryon cooling that c haracterizes the simulation s ana- 
lyzed in the work by iFang . Humphrev fc Buotc' ('200^ 
(exte nsively described in Kravt sov. Nagai & VikhlininI 

l200i 
respect to 



20051: iKravtsov. Vikhlinin fc Nagail 

Nagai. Vikhlinin fc Kravtsov! l2007l l with 
those discussed here. The stellar fraction /*(< r) = 
M*(< r)/M{< r) in the central part of our Set 1 
clus ters (i.e. < -R500) is esti m ated to be smaller than 
in iFang. Humphrevfc Buotd (|2009l ) simulations, by 
about a factor of ~ 1.5. Though, it is definitely higher 
than expected from o bserv ations of real clusters (e.g. 
iLin. Mohr fc Stanford! 120031 '). The implementation of 
cooling in GADGET-2 reduces the overcooling problem, 
naturally preventing strong rotation to get established. 

In Fig. [4] we plot the variation of vtan with redshift, 
calculated in the innermost region of g51. The peak shown 
around redshift ~ 1.5 is likely to be driven by the last major 
merger occurring to g51, and is not related to a quiescent 
build-up of mass and therefore of rotation. Instead, within 
the redshift range ~ 1.5 — 0.5, a general, although not 
smooth, increasing trend of vtan can be seen in the plot, 
which is likely to be explained as the result of the collapsing 
process under angular momentum conservation, although it 
is difficult to show it quantitatively. At lower redshift, it is 
worth to point out an interesting feature, that is the sudden 
drop of Vtan, steeply decreasing twice at 2 ~ and ~ 0.3. 
The breaks in this expected general trend are not directly 
related to any major event, and a deeper investigation of 
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Figure 4. Evolution with redshift of the tangential component 
of the ICM velocity in the innermost region (< O.liJsoo) of g51. 



the ICM internal dynamics has then been performed in 
order to understand the possible origin of this unexpected 
behavior. 

ICM velocity maps. The panels in Fig. O show the 
two-dimensional velocity field in the best equatorial plane 
in the central slices of g51. Each velocity vector has a length 
proportional to the absolute value of the velocity in that 
point of the plane. The dashed circles mark the innermost 
region enclosed within 0.1i?5oo (smaller circle) and i?5oo 
(larger circle). 

The velocity maps catch one of the two major decreases 
in the curve of vtan, in particular the one at roughly z ~ 0.3, 
which is the first significant break in the increasing trend 
shown up to redshift ~ 0.5. Clearly, one can see the passage 
of a gas-rich subhalo (thicker circle) through the best 
equatorial plane, onto which the gas velocity field has been 
projected in the Figure. The subhalo is the only gas-rich 
subhalo approaching the central region of the simulated 
cluster. 

The two panels in Fig. [5] refer to redshift z ~ 0.314 
(top) and z ~ 0.297 (bottom), and show the best moment 
right before and after the first passage of the substructure 
through the equatorial plane. The steep decrease of vtan 
does not start at this moment nor does it reach the lowest 
value, but these two redshift snapshots have been judged 
to best show a plausible explanation of the suppression of 
rotation while it is happening. In fact, from the velocity 
fields we notice that the gas shows a rotational motion with 
velocities of ~ 380km/s in the innermost region, close to 
the smaller dashed circle, while the subhalo is approaching 
(upper panel). This rotational pattern is evidently disturbed 
in the lower panel, where the subhalo has already passed 
through the plane, its gas particles get probably stripped by 
the main halo gas and contribute to decrease the velocity 
values to ~ 275km/s. 

Let us stress that there are several DM-only substruc- 
tures permanently moving within the cluster and close to 
the innermost region, but they do not disturb the build-up 
of rotation as gas-rich subhaloes do. 

The decrease of vtan at redshift ^ shows an analogous 
behavior. 
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Figure 5. Gas velocity fields at z ~ 0.314 {upper panel) and 
z ~ 0.298 (lower panel) projected onto the plane perpendicular to 
the direction of the gas mean angular momentum in the innermost 
region. The smaller and larger dashed circles mark respectively 
the regions of O.liJsoo and -R500, while the grey ones are DM- 
only subhaloes and the black circle is the gas-rich halo passing 
through the equatorial plane. The coordinates in the graphs are 
in comoving units. 



4 ROTATIONAL CONTRIBUTION TO TOTAL 
MASS 

In this section we compare the contribution coming from ro- 
tational motions that should be considered in the estimation 
of total mass with the total mass calculated for the simu- 
lated clusters. Formally, the total cluster mass M enclosed 
within a closed surface S is given by Gauss's Law 

1 

4^ , 

where $ is the cluster gravitational potential and G is the 
gravitational constant. Under the assumptions that the ICM 
is a steady-state, inviscid, collisional fluid, we can replace 
the term V<I> with the terms involving gas pressure and ve- 
locity using Euler's equation as follows: 

1 

4^ 



M 



V$ ■ d S, 



(1) 



M 



1 



Ps 



VP, - (v ■ V)v 



(2) 



where pg and Pg are the gas density and pressure respec- 
tively. While the pressure term within the integral represents 
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Figure 6. Mass profiles of the relaxed cluster g51 {upper left panel) and of the disturbed system gl (upper right panel), at z = 0. 
Top: in each panel, the radial profiles of the true mass, Mtrue (solid line), and of the component estimated from the gas rotation, Mrot 
(dotted line), out to -R500 are shown. Bottom: ratio of the two mass terms Mrot/Mtrue as function of radius for g51 (lower left panel) 
and gl(lower right panel). 



the contribution of the gas random motions, both thermal 
and turbulent, the velocity term includes the ordered mo- 
tions in the ICM, i.e. rotational and streaming motions. For 
the purpose of our work, we are mainly interested in the 
contribution to the pressure support given by the rotational 
motions of the hot intracluster gas, and we therefore sepa- 
rate the velocity term in Eq. [2] into 

and 

by evaluating (v • V)v i n spherical coordinates 
ijBinnev fc Tremain3 l2008l : iFang. Humphrey fc Buotd 
I2OO9I ). In particular, the term due to streaming motion, 
Mstr, is likely to be less relevant than Mrot, especially 
for relaxed clusters. Therefore we explicitly calculate 
the rotation term for our clusters g51 and gl at several 
redshifts. In our analysis we compare this term with the 
total "true" mass of gas and dark matter for the cluster, 
Mtrue, computed directly by summing up all the particle 
masses. 

In Fig. [6] we show the radial profiles of both Alrot and 
Mtrue, and their ratio, for g51 (left panel) and gl (right 
panel) at 2: = 0. In order to evaluate Eq. [3] we consider 
the full three-dimensional structure of the velocity field 
without any assumption of spherical symmetry, as stated 
in Gauss's theorem, so that our calculation is completely 
independent of the cluster geometry. In our approximation 
of the integral that appears in Eq. [3l the closed surface S 
has been replaced with a spherical shell at each radial bin, 
and we sum over the single particle contributions in the shell 
instead of dividing the surface in cells. Specifically, at each 
radial bin we associate an effective area to each gas particle 
in the shell and compute the integrand term considering 
the velocity (vg + v'^) of the particle. We make use of a 
radial binning up to -R500, each bin of 14/i~^kpc, consistent 



with the profiles of vtan{r) discussed in Section [3.21 

The profiles shown in Fig. [6] have similar trends for 
both g51 and gl, but one can clearly notice that the 
more significant rotational pattern found in the innermost 
region of the disturbed cluster gl with respect to g51 is 
refiected here in a more significant contribution to the total 
mass of Mrot- In general, the rotation appears to be more 
important in the innermost regions than in the cluster 
outskirts. Indeed, out to 7?5oo the rotational component of 
the total mass accounts for ~ 2% of the true mass, Mtrue, 
in g51 and slightly more in the case of gl (4%). In contrast, 
Mrot plays a more significant role at radii < O.I-R500, where 
its contribution can be up to 9 — 17% of Mtrue, the highest 
value reached in the core region of gl. Although we consider 
purely rotational motions, other non-thermal motions 
should be accounted for as well and we can conclude that, 
if rotation establishes, it can significantly contribute to the 
total pressure support to the cluster weight. 



5 EXTENDING THE STATISTICS TO A 
LARGER SIMULATED SAMPLE 

In order to gain a more reliable overview of the phenomenon 
of the appearance of rotational patterns in the innermost 
regions of cluster-like haloes, our analysis has also been 
performed for Set 2, a larger sample of simulated objects 
spanning a range of M2oom from ~ 5 x lO^*ft~^M0 to 
- 2.2 X 10^^fe"^MQ. 

The selection of the sample has been carried out such 
that we isolate all the haloes in the cosmological box with 
virial mass above a mass threshold, chosen to keep a sta- 
tistically reliable number of objects throughout the red- 
shift range explored. At 2 = 0, we selected 26 haloes with 
M2oom > 5 X lO^'*fe~^M0, for which a visual representation 
is given in Fig. [T] At higher redshift, the mass threshold is 
lower in order to have a fair sample to investigate statisti- 
cally. We calculate the distribution of the rotational veloc- 
ity on the best equatorial plane in the innermost region (i.e. 
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< O.I-R500) of each selected halo, sampling the redshift range 
[0,0.5]. 

5.1 Distribution of rotational velocities at various 
redshifts 

In Fig. [S] we show the distribution of vtan , calculated < 
O.I-R500 in the same way as for g51, for the sample of cluster- 
like haloes belonging to the Set 2. The histograms show 
how significant rotation is over the range of masses and red- 
shits, confirming the intermittent nature of the phenomenon. 
Starting form the upper-left panel to the bottom-right one, 
redshift increases from to ~ 0.5 and the sample con- 
sists of a number of objects varying between 26 and 44. 
The mass threshold chosen, Mth, is 5 x 10'^'', 4 x 10^'',3 x 
10^'', 2 X 10^'^h~^MQ for the four redshift values considered 
{z — 0,0.1,0.3,0.5, respectively). We notice that in gen- 
eral the velocity distributions are mainly centered around 
values of 200 — 300km/s, with a mild, though not substan- 
tial, shift towards higher values for intermediate redshift. 
The shaded histograms on top of the ones plotted for red- 
shifts z = 0.1,0.3,0.5 refer to subsamples of haloes with 
M2oom > 5 X lO^''/i~^M0 (the same as the one used at 
z — 0). This comparison is meant to show a clear evidence 
that also the higher-mass subsample actually agrees with 
the general trend. Although the number of objects in these 
subsamples decreases for increasing redshift, this has been 
done in order to confirm the idea that indeed the distribu- 
tion of the rotational velocities is peaked around quite low 
values. Except for some outliers, we can generally exclude 
any clear monotonic increase of the typical value of Vtan in 
the cluster innermost region (always kept to be < O.I-R500) 
associated to the assembling of the cluster-size haloes. In 
our simulations we cannot find any quiescent build-up of 
rotation as a consequence of mass assembly, and the distri- 
bution among the volume-selected sample is not dramati- 
cally changing with redshift. 



6 DISCUSSION AND CONCLUSION 

In this work, we have presented the result of a study over 
two sets of hydrodynamical simulations performed with the 
TreePM/SPH code GADGET-2. The simulations include 
radiative cooling, star formation, and supernova feedback 
and assume slightly different cosmological models, a ACDM 
one and a WMAP3 one. The main target of this analysis has 
been the importance of rotational gas motions in the central 
regions of simulated cluster-like haloes, as it is thought to 
be a crucial issue while weighing galaxy clusters and iden- 
tifying them as relaxed systems. The objects selected from 
our samples guarantee a wide variety of virial masses and 
dynamical structures, so that a reliable investigation of this 
phenomenon is allowed. 

Our main results can be summarized as follows: 

• As main conclusion, we notice that the occurrence of 
rotational patterns in the simulated ICM is strictly related 
to the internal dynamics of gas-rich substructures in a com- 
plicated way, so that it is definitely important to take it into 
account as contribution to the pressure support, but it's not 
directly nor simply connected to the global dynamical state 
of the halo. 



• In the first part of our analysis we focused on g51, a sim- 
ulated cluster with a very smooth late accretion history, iso- 
lated and characterized by few substructures in comparison 
to the other massive objects within Set 1. Also, we compare 
it with a highly disturbed system (gl). Even in the radiative 
simulation of this cluster, likely to be considered relaxed in 
a global sense, no clear rotation shows up at low redshift be- 
cause of some minor merging events occurring close to the 
innermost region: the rotation of the core is found to be 
an intermittent phenomenon that can be easily destroyed 
by the passage of gas-rich subhaloes through the equatorial 
plane. Gas particles stripped from the subhalo passing close 
to the main-halo innermost region (< O.I-R500), are likely to 
get mixed to the gas already settled and contribute over few 
orbits to change the inclination of the best equatorial plane, 
suppressing any pre-existing rotational pattern. 

• The velocity maps plotted in Fig. [S] show several DM- 
only subhaloes moving close to g51 central core. In our study, 
they have been found not to disturb in any significant way 
the ordered rotational gas motions created in the innermost 
region. The central gas sloshing is mainly set off by gas-rich 
subhaloes, especially if they retain their gas during the early 
passages through the core. Interesting work on numerical 
simulations have been found to b e relevant for the result 
presented here, as the study from lAscasibar fc MarkevitchI 
( 2006) on the origin of cold fronts and core sloshing in galaxy 
clusters. 

• Mass measurements based on HEH are likely to mises- 
timate the total mass of galaxy clusters because of contri- 
butions by non-thermal gas motions that have to be consid- 
ered. In agreement with previous works, we also find that 
significant rotation of the ICM can contribute to the pres- 
sure support. While several studies have been carried out on 
turbulent motions in the ICM and on th eir effect on the clus- 
ter mass estimates (e.g. [ Rasia. Tormen fc Moscardini 2004 
Fang. Humphrev fc Buot' 




Lau. Kravtsov fc Nagai 



ol 



20091 : IZhurayleva et al.1 120101), ori ly lately 
work by iFang. Humphrey fc Buotd (|2009l ') 
iLau, Kravtsov fc Nagail (|2009 i) have been ad( 



the 
and 
addressing 

the ordered rotational patterns that could establish in 
the innermost region ICM as the result of the cluster 
collapse. Therefore, a comprehensive analysis of the de- 
tails of rotation build-up and suppression both in single 
high-resolution case-studies and in larger, statistically 
significant samples is extremely interesting, especially for 
relaxed objects where this should be more important than 
turbulence. Focusing on rotation specifically, we calculated 
the corresponding mass term, Mrot, for the two clusters g51 
and gl. As expected from the tangential velocity profiles at 
redshift z = 0, the mass term coming from ICM rotational 
motions contributes more in the case of gl than in g51, 
providing evidence that rotational support of gas in the 
innermost region is more significant in the former than in 
the latter. While Mrot accounts for few percents at radii 
close to -R500 in both cases, in the central regions up to 
~ 17% of the total true mass in gl is due to rotational 
motions of the ICM. As regards g51, this contribution is 
less important, as no strong rotation has been found at 
z = 0, but it still reaches a value of ^ 10% for the pressure 
support in the cluster core. 

• Extending the analysis to a larger sample, we have in- 
vestigated the statistical distribution of rotational velocity 
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Figure 8. Distribution of the rotational velocity for a sample of simulated clusters extracted from the Set 2. The tangential component 
of the gas velocity, vtan, is the mean value within O.I-R500, calculated in the plane perpendicular to the mean gas angular momentum 
in the same region. Top left to bottom right, the panels refer to a sample of objects from the Set 2 at different redshifts, sampling the 
range [0,0.5]. The sample of objects has been selected according to a threshold in mass, M^^, which is lower at higher redshift in order 
to keep a fair statistics. While M^jj is different at each redshift, the shaded regions represent subsamples of haloes selected to be more 
massive than 5 X 10^'*/i~^Mq, as a comparison to the distribution at 2 = 0. 



over dynamically-different clusters, isolated in a limited- 
volume simulated box such that their virial mass (A'/2oom) 
is above a chosen threshold. At 2 = as well as at higher 
redshift up to ~ 0.5, a fair sample of cluster-size haloes let us 
infer that, on average, no high-velocity rotational patterns 
show up in the halo cores (i.e. in the region < O.l-Rsoo)- 
Also for the clusters of Set 2, we find typical values of 
~ 200 — 3Q0km/s for the rotational velocity in the inner- 
most region. 

• We do not find any increasing trend of the rotational 
velocity distribution peak with decreasing redshift, that can 
correspond to the smooth mass assembly of the cluster-like 
halo through collapse. Although such trend is generally ex- 
pected, it must be easily suppressed by internal minor events 
disturbing the halo central region. 

We conclude that the build-up of rotational patterns in the 
innermost region of galaxy clusters is mainly related to the 
physical processes included in the csf run to describe the 
intracluster gas. On the contrary, numerical e ffects such as 
differ ent implementations of artificial viscosity (jPolag et al.l 
[iooi) do not affect in any significant way our results (see 



Appendix |X1 for a detailed discussion). 

An analogous conclusion can be drawn with respect to 
the differences between the two samples introduced by cos- 
mology and resolution. For both Set 1 and Set 2 the build- 
up and suppression of rotational patterns in the halo cen- 
tral part is found to be mainly related to the physics in- 
cluded in the radiative run. In fact, comparable subsamples 
of the two sets in the csf simulations show very similar dis- 
tributions of rotational velocities for the ICM component in 
the halo innermost region, meaning that the shape of the 
distribution is essentially dominated by the physics of the 
gas. Usually, relaxed clusters are assumed to have lit- 
tle gas motions. Therefore they are likely to be the best 
candidates for the validity of the HEH, on which mass esti- 
mations are based. Nevertheless, rotational motions should 
establish preferentially in relaxed clusters with respect to 
disturbed systems as a consequence of the assembling pro- 
cess, potentially representing a danger for relaxed cluster 
masses. Here, however, we find that the processes described 
in the paper save the reliability of the HEH-based mass 
determinations in most of the cases. In fact, rotational mo- 
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tions are not significant enough to compromise dramatically 
mass determinations with the exception of few outliers. In 
our simulation, the identification of relaxed or non-relaxed 
clusters according to the presence of gas rotation in the cen- 
tral region is not straightforward, since it has been shown 
to appear and disappear periodically. Its contribution has 
to be considered whenever is present, but it is not directly 
related to the global state of the simulated halo. Also, it is 
likely to be strongly influenced by the overcooling problem 
affecting hydrodynamical simulations, which has the effect 
to enhance the process of building up rotational patterns in 
the ICM in the innermost regions of simulated clusters. 

Although various theoretical and numerical studies 
in addition to the present work have been investigat- 
ing the existence of gas bulk, non-thermal motions and 
the possible ways to detect them in galajcy clusters (e.g. 
Fang. Humph rey &LB uotc 200J^; iLau. Kravtsov fc Nagail 
20091 : iZhuravleva et al.ll2010f). little is known from o b serva - 



tions. In a recent study, Lagana. de Souza fc Kellerl (|2009l ) 
have made use of assumptions from theoretical models and 
numerical simulations about cosmic rays, turbulence and 
magnetic pressure to consider these non-thermal contribu- 
tions to the total mass measurement for five Abell clusters. 
From a pure observational point of view, previous work has 
been able to confirm only indirect indications of bulk gas 
motions associated to merging events in galaxy clusters (see 
iMarkevitch fc Vikhlinin||2007 , for a review) or evidences for 
turbu lent gas motions, likethe ones found in the Coma clus- 
ter in lSchuecker et al] (|2004l ) or those inferred, on the scale 
of smaller-mass systems, from the effects of resonant scat- 
tering in the X-ray emitting gas eous haloes of large ellipti- 
cal galaxies (| Werner et al .1120091 ). Although not possible so 
far, the most direct way to measure gas motions in galaxy 
clusters would be via the broadening of the line profile of 
heavy ions (like the iron line at --^ 6.7keV in X-rays) for 
which the expected linewidth due to impact of gas mo- 
tion is much larger than the width due to pure thermal 
broadening. The possibility to use the shape of the emis- 
sion lines as a source of information on the ICM velocity 
field a s bee n discussed in detail in Inogamov fc SunvaevI 
l|2003l) and ISunvaev. Norman fc BrvanI and lately 



Rebusco et al.l ( 20081 ). Though, the investigation of the 



imprint of ICM motions on the iron line profile requires 
high-resolution spectroscopy, which will become possible in 
the near future with the next-generation X-ray instruments 
such as ASTRO-H and IXO. This will aUow us to directly 
detect non-thermal contributions to the cluster pressure 
support, such as rotational patterns in the ICM, and en- 
able us to take this correctly into account as contribution 
to the total mass estimate. Ultimately, this is likely to be 
an important issue to handle in order to better understand 
deviations from the HEH, on which scaling laws are usually 
based. 
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APPENDIX A: 
VISCOSITY 



EFFECTS OF ARTIFICIAL 



The runs studied in the present work are the csf simula- 
tion, including radiative cooling, star formation and super- 
nova feedback, and the non-radiative simulation (labelled 
as ovisc), where the original parametr ization of artificial vis- 
cosity bv lMonaghan fc Gingoldl ()l983[ ) is used. We comment 
here on the effects of the artificial viscosity on our study by 
considering two additional non-radiative runs of our simu- 
lations, carried out with alternative implementations of the 
artificial viscosity scheme. In particular, we label as svisc 
the non-radiative run with slightly less numerical viscosity 
based on the signal velocity approach of Monaghan (1997), 
and as Ivisc the modified artificial viscosity scheme where 
each particle evolves with its own time -dependent viscosity 
pararn eter (as originally suggested in iMorris fc Monaehanl 
l|l997ll ). For a detailed descri ption of these non-radiative 
runs we refer the reader to Dolag et al.l ([2005). 

In Fig. lAll the distribution of the value of vtan for the 
Set 1 is presented for the different runs. With the solid 
black line we refer to the csf simulation, while the other 
histograms represent the non-radiative runs: ovisc (dotted 
red line), svisc (green dashed) and Ivisc (blue, dot-dashed). 
As one can see from the Figure, the csf simulation shows a 
different distribution of rotational velocity in the innermost 
regions of the cluster-like haloes with respect to the non- 
radiative runs, in which the difference in the implementation 
for the artificial viscosity does not produce significant differ- 
ences in the three distributions. In Fig. lAll is evident that 
the difference between the radiative and non-radiative runs 
is much larger than the difference among the non-radiative 
runs themselves. All the non-radiative runs similarly show 
that the largest fraction of clusters have very low values of 
Vtan, and differ mostly in the lack of a high- velocity end of 
the distribution from the csf case. The cooling of the core 
allowed by the physics included in the csf run is plausi- 
bly responsible for the presence of a significant fraction of 
clusters with high velocity values, which do not exist in the 
non-radiative runs. 

From Fig. |XT] we can confirm that, accordingly to what 
is expected, the overcooling problem coming along with nu- 
merical simulations of galaxy clusters leads to a more signif- 
icant build-up of rotation in the core, refiected in an overall 
shift towards higher values of the distribution of Vtan - 

For the purpose of our study, we can conclude that the 
main effects in the establishment of rotational patterns in 
the central region of simulated clusters are introduced by 
the physical processes describing the gaseous component, 
included in the csf simulation. The amount of turbulence, 
referring in particular to small scale chaotic motions, was 
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Figure Al. Distribution of rotational velocity (calculated in the 
region < O.lHsoo) for the Set 1 in the redshift range [0 — 0.5]. 
The different runs are plotted: csf {solid, black line with trian- 
gles), ovisc (dotted, red with crosses), svisc (green, dashed with 
diamonds), Ivisc (blue, dot-dashed with asterisks). 
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tigated here (e.z. iDolag et al.]|2005f ). However, the effect of 
artificial viscosity on rotation does not depend strongly on 
the specific numerical scheme implemented, and we safely 
confirm our conclusions about the more significant effect 
that minor events occurring close to the cluster central re- 
gion have on the survival of gas rotational motions. 



APPENDIX B: ELLIPTICITY OF ICM 



As in terestingly suggested in a recent work by iLau et al] 
(|2010l ). certain observable features of the cluster X-ray emis- 
sion, like a fiattening of cluster shapes, could unveil the pres- 
ence of rotationally supported gas, especially in the inner- 
most region. Referring to the X-ray surface brightness maps 
in Fig. [TJ we compare here the ellipticity profiles extracted 
from the maps in the three projections to the rotational ve- 
locity profile, for the two cases of study extracted from Set 
1 (namely, g51 and gl) at redshift z — 0. The ellipticity was 
calculated by means of the IRAF task ellipse and the re- 
ported error bars are those obtained by the fitting method. 
In Fig. IBll the vtan profile is marked by the solid curve for 
each cluster, while the dotted, dashed and dot-dashed lines 
refer to the ellipticities from the three projected maps. Let 
us note that the ellipticity profiles do not extend to the 
very central part, since the complex structure in the inner- 
most region does not allow for a simple determination of 
ellipticities. We can conclude that only a mild relation be- 
tween ellipticity and rotation is found in these clusters and 
the difference in the trends shown by the rotational velocity 
profiles is stronger than the difference among the elliptic- 
ity profiles of the two clusters. Nevertheless, comparing the 
csf run (upper panels) to the ovisc run (lower panels) we 
can confirm that the cluster shape is actually marked by 
the physics included in the radiative run, especially in the 
gl case, where some rotation establishes in the csf simula- 
tion. In conclusion, gas cooling can affect the gas shape by 
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Figure Bl. Radial profiles of ellipticity {left y-axis; dotted, dashed and dott-dashed curves) for the three projection axis compared with 
the radial profile of the gas rotational velocity {right y-axis; solid curve). The panels refer to the csf and the ovisc simulations {upper 
and lower panels, respectively) of g51 {left) and gl {right) at z = 0. 



flattening the X-ray isophotcs but wc do not find a strong 
significant evidence for that in our case-study clusters, in 
agreement with the milder rotation found and the nature 
of such rotational patterns (Uke in the gl example), which 
are likely to be temporary effects driven by major merging 
events. 
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